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We classify all spherically symmetric and homothetic spacetimes that are allowed kmematically 
by constructing them from a small number of building blocks. We then restrict attention to a 
particular dynamics, namely perfect fluid matter with the scale-free barotropic equation of state 
p = an where < a < 1 is a constant. We assign conformal diagrams to all solutions in the 
complete classification of Carr and Coley, and so establish which of the kinematic possibilities are 
realized for these dynamics. We pay particular attention to those solutions which arise as critical 
solutions during gravitational collapse. 
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I. INTRODUCTION 



Self-similar and spherically symmetric solutions play an 
important role in many areas of physics. In Newtonian 



physics, using suitable space and time coordinates x and 
t, self-similarity is an ansatz f{x,t) — folx/t) for some 
field / that is of interest. On the one hand, such an 
ansatz simplifies the field equations by eliminating one 
coordinate. (In spherical symmetry, in particular, the 
self-similar ansatz simplifies field equations that are par- 
tial differential equations in r and t to ordinary differen- 
tial equations in one variable r/t.) On the other hand, 
self-similar solutions are often attractors, or intermedi- 
ate attractors, and are therefore good approximations to 
generic solutions during some stage of their evolution. 

A classic example is the solution of the Riemann problem 
in 1-dimensional fluid dynamics, where at i = the fluid 
has one constant velocity, density and entropy for x < 
0, and another set of constant values for a: > 0. The 
resulting shock solution is self-similar in the coordinate 
x/t. The reason for self-similarity is that the initial data 
do not contain any length scale. The self-similar solution 
is a good approximation over an intermediate range of 
scales, larger than the molecular scale where the fluid 
assumption breaks down but smaller than the scale on 
which the solution varies smoothly. 

In general relativity, the most common definition of self- 
similarity is homotheticity: the spacetime admits a vec- 
tor field with the property that the Lie-derivative of 
the metric along this field is a constant times the metric: 



(1) 



(The choice of the constant as —2 is a convention, the 
negative sign being chosen so that points towards in- 
creasing curvature.) This definition expresses the scale- 
invariance of the spacetime metric, and is analogous to 
the Newtonian notion of self-similarity. With spherical 
symmetry, it is possible to introduce coordinates t and 
r so that all metric coefficients depend only on r/t. In 
general relativity, self-similarity has an important appli- 
cation in gravitational collapse: during collapse whose 
initial data are close to the black hole threshold, certain 
self-similar "critical" solutions act as universal interme- 
diate attractors 
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The continuous self-similarity (CSS) that we have dis- 
cussed already can be generalized to a discrete self- 
similarity (DSS), in which f{x,t) depends arbitrarily on 
x/t but can also depend on ln|i| in a periodic manner, 
with some period A. Clearly such solutions are scale- 
periodic rather than scale-invariant. Several known criti- 
cal solutions in gravitational collapse, in particular those 
with scalar field matter and in vacuum gravity, are DSS 
rather than CSS. In this paper we are interested in per- 
fect fluids, where the critical solutions are CSS, and from 
now on we restrict our discussion to this case. We refer 
to CSS when we specifically mean the continuous sym- 
metry, and to self-similarity in any statements that are 
also true for DSS. 

We proceed in three steps. We first discuss the kinematic 
consequences of CSS in spherical symmetry, generalizing 
the discussion in Ori and Piran II and Nolan (|. We 
then discuss the dynamics with perfect fluid matter and 
examine which of the kinematic possibilities are realized 
dynamically. We assume the linear barotropic equation 
of state p — a/i, which is the only one that is compatible 
with exact self-similarity. We use the classification by 
Carr and Coley |^ of all spherically symmetric CSS per- 
fect fluid solutions to draw the corresponding conformal 
diagrams. Finally, we look more closely at those solutions 
that appear as critical solutions in gravitational collapse. 



II. KINEMATICS OF SPHERICALLY 
SYMMETRIC AND HOMOTHETIC 
SPACETIMES 



In this section we discuss the kinematics of spherical sym- 
metry, without using the Einstein equations. Any spher- 
ically symmetric spacetime is naturally the product of a 
two-dimensional reduced spacetime crossed with the 
2-sphere S*^. Any point in the reduced spacetime corre- 
sponds to a 2-sphere of area AttR^ . Here R transforms as 
a scalar on the reduced spacetime, and R — is a bound- 
ary of the reduced spacetime, which can either be a regu- 
lar center or a central curvature singularity. A local mass 
function m can be defined by 1 — {2m/R) = V^/iV^i?. 
Any point in the reduced spacetime where V^i? is time- 
like is a closed trapped surface in the full spacetime and a 
surface where V^i? is null is called an apparent horizon. 

We now restrict attention to spherically symmetric CSS 
spacetimes. The integral curves of the homothetic vec- 
tor field (homothetic lines) provide a natural fibration 
of the 4-dimensional spacetime, and therefore a natural 
foliation of the reduced spacetime. Labelling the inte- 
gral curves provides a natural global coordinate on the 
reduced spacetime, which we shall call z. In this paper 
we construct conformal diagrams for the reduced space- 
time. Apart from any singularities and conformal bound- 
aries, lines of geometric significance include the homoth- 
etic lines, the lines of constant R and the fluid world 



lines. 

We shall describe two building blocks from which any 
spherically symmetric CSS spacetime can be built. Our 
arguments in this section are purely kinematic. They are 
therefore valid for any matter. On the other hand, for 
a given type of matter not every spacetime that can be 
assembled from the building blocks can be realized as 
a solution of the field equations. We illustrate this in 
Section IV for perfect fluid matter. 



The key kinematic consequence of self-similarity in gen- 
eral relativity is the existence of a strong curvature sin- 
gularity. It can be shown from (|l|) that in the direction 
^'^ any homothetic line runs into a curvature singularity 
in finite proper time, proper distance or affine parameter 
distance s, unless the curvature vanishes identically on 
the curve [|[ |^. If the origin of s is chosen so that the 
singularity is at s = 0, the Kretschmann scalar scales as 
s~*, the Ricci scalar as and the mass as s~^. Ev- 
ery homothetic line is infinitely extended in the direction 
— (i.e. with increasing s). 

Any homothetic spacetime has therefore the manifold 
structure of the half-line < s < oo crossed with an 
interval in z (which can be finite or infinite), s = is 
always a curvature singularity and s = cxd is always a 
physical infinity. The boundaries z = Zmin and z = 2,nax 
can be either regular centers or central singularities. The 
manifold structure is shown in Fig. |l|. It is important to 
note that this on its own does not determine the causal 
structure, but it does mean that all components of the 
singularity are connected. The usual conformal diagram 
of a closed Friedmann universe with disjoint big bang 
and big crunch singularities, for example, cannot be re- 
alized in a spherically symmetric CSS spacetime. Never- 
theless, fluid world lines can begin on one section of the 
singularity (big bang) and end on another (black hole/big 
crunch). The difference is that in a spherically symmet- 
ric CSS solution these two components of the singularity 
must be connected either directly or through a null sin- 
gularity. 

We obtain our two building blocks from considering ra- 
dial null geodesies that are also homothetic lines. Ori 
and Piran ^ call such a line a simple radial null geodesic 
(SRNG). Assuming that the spacetime contains an iso- 
lated SRNG, with neighboring homothetic lines are time- 
like on one side of the SRNG and spacelike on the other, 
there are two possibilities: either the homothetic lines 
converge with respect to radial null geodesies as they ap- 
proach the singularity at s = or they diverge. 

The first possibility is illustrated in Fig. ^. We call this 
building block a "fan" . The homothetic lines converge as 
s —^ 0. The singularity s = 0, for an interval of z, is then 
just a single point. Any non-simple radial null geodesic 
is given by z = z{s). It is clear from the figure that z(s) 
is repelled from zq on either side as s — > 0. Conversely, 
z{s) is attracted to zg as s ^ oo. (z = zo,s = oo) 
is therefore not a point but an extended piece of null 
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s=0 s=0 

FIG. 1: The manifold structure of any spherically symmetric 
self-similar spacetime: an interval Zmin < z < Zmayi times the 
half-line < s < oo. Note that this diagram is not a confor- 
mal diagram: although s = is shown here as a horizontal 
line, it could be a single point or a line with spacelike, time- 
like or null segments. Homothetic lines z = const, are shown 
bold with arrows that point to the singularity s = 0. The 
four cases where none, one or both of the limits z = Zmin and 
Zmax are singularities are shown. 



infinity. In summary, the spacetime contains a fan when 
it contains a SRNG z = zq such that z(s) for non-simple 
radial null geodesies approaches zq as s — > oo, that is 
it obeys [z{s) — zq] dz/ds < 0. A fan implies a point 
singularity at s = for a finite range of z, and a null 
infinity at (s ~ oo, z = zq). 

The opposite structure is illustrated in Fig. ||. Here the 
non-simple radial null geodesies z = z{s) approach z — zq 
as s ^ 0, and are repelled from zq as s — > oo. Therefore 
the singularity (z = zq, s = 0) is reached not just by one 
but many radial null geodesies, and so it is a null line. 
Conversely, s = oo for a range of z on each side of zq is a 
just a single point, rather than an extended piece, of null 
infinity. We call this structure a "splash" because the 
homothetic line z = zq runs into the singularity at s = 
and then bifurcates to form the null singularity. The 
spacetime contains a splash when it contains a SRNG 
z = Zq such that z(s) for non-simple radial null geodesies 
approaches zq as s ^ 0. A splash implies a null central 
singularity at s — Q. 

Note that the coordinates s and z break down at both 
s — Q and s = oo in a fan or splash. Regular double 
null coordinates for a fan and a splash are constructed 
in Q . (This construction also proves that the singularity 
s = in a splash is null, a fact which was asserted but 
not proved in jl] and |^.) Note also that in principle the 
homothetic lines could be timelike (or spacelike) on both 
sides of the SNRG. In this case, one would only have the 
upper (or lower) half of the diagrams shown in Fig. |^ and 
Fig. ^, corresponding to what one might term a "semi- 
fan" or "semi-splash" respectively. However, we do not 
show the corresponding figures explicitly. 




singularity s=0 



FIG. 2: A fan: a range of homothetic lines meet in a point 
on the singularity. Homothetic lines are shown bold with 
arrows that point to the singularity, in the direction One 
of them, z = zo, IS also a radial null geodesic (a SRNG). 
Four other (non-simple) radial null geodesies z — z{s) are 
shown as thin lines. The conformal diagram on the left is not 
compactified. On the right, the same conformal diagram has 
been compactified to show that {z — zo,s = oo) is a piece of 
null infinity. 




FIG. 3: A splash: A SRNG z ^ zo is flanked by homothetic 
lines that diverge from it as the singularity is approached. 
Homothetic lines z = const, are shown bold with arrows that 
point to the singularity. The SNRG branches and runs along 
the singularity. The singularity s = is null and is at finite 
distance. The conformal diagram on the left is not compact- 
ified. The conformal diagram on the right is compactified to 
show that a range of homothetic lines meet in a point at null 
infinity. 

III. PERFECT FLUID MATTER 

We now examine spherically symmetric CSS spacetimes 
with perfect fluid matter with the equation of state 
p — afi, where p is the pressure, p, is the total energy 
density, and < a < 1 is a constant. We briefly intro- 
duce notation based on ||^ and The metric is given 
in comoving coordinates by 

ds2 ^ -e^'^dt^ + e^^dr^ + i?2 (2) 

where < r < oo labels fluid worldlines. The fluid ve- 
locity is therefore = (e"", 0, 0, 0). Spherically sym- 
metric CSS solutions can be put into a form in which all 
dimensionless quantities are functions only of the dimen- 
sionless self-similar variable z = r/t. In particular, we 
enforce homotheticity by demanding that 

v{t,r) = v{z), X(t,r)=X{z), 
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R{t, r) = rS{z), 8TTGfi{t, r) = r'^-qiz). (3) 
The homothetic vector defined by (|^) is 



^„ d d d 

^ d'' dt dr 



(4) 



in these coordinates. Varying z for a given t specifies the 
spatial profile of various quantities. For a given value of 
r (i.e., for a given fluid element), it specifies their time 
evolution. The remaining gauge-freedom in CSS solu- 
tions corresponds to a scaling of the r and t coordinates. 

Two quantities have special physical significance ^ . The 
first is the CSS mass function, which is defined in terms 
of the Hawking mass m{r,t) and the area radius R{t): 



Miz)^^. 



(5) 



One has an apparent horizon wherever M — 1/2. The 
second is the CSS velocity function 



V{z) 



e^-'^z. 



(6) 



which represents the velocity of the spheres of constant 
z relative to the fluid. In geometric terms. 



1 



(7) 



This should not be confused with the velocity of the fluid 
with respect to the lines of constant R, the "radial 3- 
velocity" , which we denote by Vr. Special significance is 
attached to values of z for which \V\ — y/a and \V\ = 
1. The first corresponds to a sonic point, the second 
to a black-hole event horizon or a cosmological particle 
horizon. A useful relation is that for radial null geodesies 



dlnr 



din 2 
l±V{z)' 



(8) 



where the + and - signs correspond to ingoing and out- 
going geodesies, respectively. This again shows that lines 
z = zo where ^^(^o) — ±1 are both radial null geodesies 
and homothetic lines, that is they are SRNGs. 

It should be noted that finite values of z sometimes cor- 
respond to zero or infinite R and zero values of z some- 
times correspond to non-zero R. Also some solutions 
span both negative and positive values of t, which means 
that z = r/t jumps from — oo to -l-oo. (This occurs at 
t — 0, which is the homothetic line orthogonal to the 
fluid worldlines.) 

Two of the Einstein equations can be solved in closed 
form to give v and A in terms of z, S{z) and ri{z). One 
then obtains a second-order ODE for S and a first-order 
ODE for T] in the independent variable z. The equation 
for cPS/dz'^ can be written as two first-order equations 
for dS/dz and dM/dz. One then obtains the system 



dM 
dz 



aiS,M,r^,z), 



(9) 



^ = b{S,M,7j,z), 



dr] 
dz 

subject to the constraint 



gjS, M, 7], z) 
V'^{S, ?7, z) — a' 



c{S,M,T],z) = 0. 



(10) 
(11) 

(12) 



Here a, 6, c, V"^ and g are known regular functions of M , 
S, 7] and the independent variable z. The constraint can 
in principle be used to eliminate z from the right-hand 
sides, so that one obtains a 3-dimensional dynamical sys- 
tem. Therefore, for a given value of a, the space of solu- 
tions is locally 2-dimensional. In practice the constrained 
and non-autonomous form is retained for calculations. 

Points where V ~ ±^/a need particular attention. In 
phase space the condition = a corresponds to a 2- 
dimensional surface, which we call the sonic surface. Eqn 
(10) implies that solutions can be continued through the 
sonic surface only on the 1-dimensional line in phase 
space where both — a and g = 0; we call this the sonic 
line. We use the term sonic point for the 3-dimensional 
surface in spacetime where = a. 

The details are summarized for completeness in Ap- 
pendix Here we discuss the effect that the existence of 
sonic points has on the local dimensionality of the solu- 
tion space. The sonic line can be divided into three parts. 
It has a "focal" segment that cannot be crossed by any 
solution curves and is a repeller. It also has a "saddle 
segment" , each point of which is crossed by two isolated 
analytic solutions. Finally, it has two "node" segments, 
each point of which is crossed by a 1-parameter fam- 
ily of solutions and an isolated solution. All the 
solutions are tangential to each other at the sonic line. 
Therefore, if one does not require the solution to be an- 
alytic, any of the solutions on one side of the sonic 
surface can be continued by any of the solutions on 
the other side. Like a saddle, a node is also crossed by 
two analytic solutions. One of these is a member of the 
family, while the other is the isolated solution. 

The global dimensionality of the families of solutions that 
we discuss below can be determined as follows. If Zmin is a 
singularity, the space of solutions locally has two free pa- 
rameters. If it is a regular center, it obeys one condition 
and so the solution locally has only one free parameter. 
If the solution crosses the sonic point analytically (i.e. at 
a node or a saddle), it loses one parameter. For example, 
critical solutions have a regular center and are analytic at 
the sonic point, so there are no free parameters, indicat- 
ing that there are at most isolated solutions with those 
properties. On the other hand, if the solution crosses a 
node as a solution, the continuation is not unique and 
so gains one parameter. The same argument applies to 
any subsequent sonic point. 

Note that the presence of a sonic point cannot be inferred 
from the conformal diagrams presented below. Each 
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sonic point is a timelike homothetic line, but a region 
in which the homothetic lines are timelike need not con- 
tain a sonic point. 



IV. CAUSAL STRUCTURE OF PERFECT 
FLUID CSS SOLUTIONS 

Carr and Coley |^ (from now on CC) classify solu- 
tions in terms of their behavior in the asymptotic limit 
\z\ — > oo, since this can only take one of a few sim- 
ple forms: asymptotically Friedmann, asymptotically 
Kantowski-Sachs, what they term asymptotically 'quasi- 
static', and asymptotically Minkowski. There is also a 
family which asymptotes to the Minkowski form at finite 
z because R ^ oo there. In the conformal diagram, all 
these asymptotic behaviors are associated with the region 
near i°. CC also discuss the behavior of solutions in the 
limit z — > and find that they are either exactly static 
or asymptotically Friedmann (corresponding to a 'regu- 
lar' origin). There are also solutions which have their 
origin at finite z because i? — > there, corresponding to 
a spacelike singularity. In the conformal diagram, these 
asymptotic behaviors are associated with the regions near 
j+ or . 

CC place considerable emphasis on the form of the scale 
factor S{z) since this specifies the physical behavior of 
the fiuid elements. They also emphasize the form of the 
velocity function V{z) since this identifies event and par- 
ticle horizons, sonic points and physical singularities. Al- 
though the V{z) representation does not yield a complete 
understanding of solutions - one is projecting solutions 
in a 3-dimensional space onto a particular 2-dimensional 
plane, so that many physically distinct solutions may be 
superposed - it uniquely specifies features relevant to the 
conformal structure. In particular, = 1 corresponds 
to a null line and infinite \V\ corresponds to a singular- 
ity. Integrating Eq. (8) shows that the value of r at 
V = 1 generally diverges for any outgoing photon, so 
this corresponds to Our strategy here is to present 
the form of the functions S{z) and V{z) for each of the 
types of solution and then to use these to construct the 
conformal diagram. We will also show some fluid lines 
and homothetic lines in each case. 

Note that we can trivially set i — )■ — i and hence z —z 
in any CSS solution, and so obtain its time reverse. As 
a matter of convention, we shall present each solution in 
only one time orientation, with "big bang" singularities in 
the past and "black hole" singularities in the future. For 
example. Fig. ^ describes an open universe that starts 
with a big bang and expands forever. Turned upside 
down, the same solution describes gravitational collapse 
from a regular initial state [Q . 

Carr et al. 0, § (from now on CCGNU) classify so- 
lutions in terms of the global causal properties of the 
homothetic vector ^'^ and this is particularly relevant 



in constructing the conformal diagram. The solutions 
are found to have five possible forms. There are type-T 
solutions for which is always timelike; these develop 
a shock, since their orbits necessarily end at an irreg- 
ular sonic point. There are type-S solutions for which 
always is spacelike; these always have \V\ > 1 and 
include some asymptotically Friedmann and asymptoti- 
cally quasi-static solutions. There are type-TS solutions, 
such as the flat Friedmann one, for which changes 
causality once. Finally there are are type-TST, STST, 
and STSTS solutions with two, three and four causality 
changes, respectively. 

Here we make this shorthand description more detailed 
by denoting the fans and splashes that separate the S and 
T regions by f and p. Finally we denote the two bound- 
aries Zmin and Zmax by either r for a regular center or s 
for a singularity. Note that the upper case letters S and 
T denote regions, while the lower case letters f, p, r and s 
denote lines. All kinematically possible CSS spacetimes 
can be constructed systematically from an alphabet of 
six letters subject to the following rules: The flrst and 
last letter are r or s. The other letters are f, p, T, S. 
Capital letters and lower case letters must alternate. 



A. Asymptotically Friedmann solutions 

There are two 1-parameter families of asymptotically 
Friedmann solutions as jzj — > oo, one with z > and the 
other with z < 0. The parameter characterizing these 
solutions measures the underdensity or overdensity rel- 
ative to the flat Friedmann solution and is also associ- 
ated with the asymptotic energy E. The z > solutions 
correspond to inhomogeneous models that start from an 
initial big bang singularity at z = oo (t = 0) and then, as 
z decreases, either expand to infinity or recoUapse. The 
z < solutions are just the time-reverse of these but we 
do not discuss them explicitly. Depending on the value 
of E, two qualitatively different types of solution can be 
distinguished: 

1) All solutions which are underdense {E > 0) or insuf- 
ficiently overdense {E < but exceeding some negative 
critical value i^crit) reach the sonic surface at V = y/a. 
They have ^ = 1 at some point zi and this corresponds 
to the cosmological particle horizon. The solutions in 
this class are of type ST for z > and are illustrated in 
Fig. ^. Those which reach the sonic line in the range of 
z corresponding to nodes may be attached to a subsonic 
{V < \/a) solution which is regular solution at the origin. 
These solutions are also described by a single parameter 
and this is a measure of the density at the origin z — 0, 
so the transonic solutions represent density fiuctuations 
that grow at the same rate as the particle horizon. They 
are generally non-analytic at the sonic point. While there 
is a continuum of regular underdense solutions, regular 
overdense solutions only occur in narrow bands (with just 
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z=0 




z=oo 

FIG. 4: The form of S{z) and V{z) for the asymptotically 
Friedmann solutions which expand forever, together with the 
conformal diagram, nomothetic lines are shown, together 
with a single fluid world line (shown broken). The arrows, 
here and in all diagrams that follow, indicate the direction of 
time. The causal structure is "singularity, region where the 
homothetic lines are spacelike, fan, timelike region, regular 
center" , or sSfTr in our shorthand notation. 



one solution per band being analytic). The overdense 
solutions exhibit oscillations in the subsonic region and 
solutions with larger numbers of oscillations form ever 
narrower bands within the one-oscillation band (in terms 
of the ranges of z at the sonic line) . The higher bands are 
all nearly static near the sonic point, although they de- 
part from the static solution as they approach the origin. 
The conformal diagram for these solutions is essentially 
identical to that for flat Friedmann itself. 

2) Solutions that are sufhcicntly overdense with respect 
to the flat Friedmann solution (i.e. for E less than the 
critical value i?crit) have the feature that V reaches a 
minimum and then rises again to infinity as z decreases 
to some finite value zs- This indicates the formation of 
a singularity at which — > and /i — > cx) for finite a 
value of z. The solutions in this class are of type STS or 
S and are illustrated in Fig. |5| and Fig. ^, respectively. 
They correspond to black holes growing at the same rate 
as the Universe. [Note that the type STS solution with 
the smallest value of zs is the one for which the minimum 
value of is 1/ y/a and this solution must touch the sonic 
surface at the value of z associated with the saddle/node 
transition.] Providing E exceeds some other critical value 
i?*, the minimum of V is below 1, as illustrated in Fig. ||, 
so there is a black-hole event horizon at Z2 and a cos- 
mological particle horizon at zi where V = 1. However, 
it should be stressed that the conformal diagram is very 
different from that for a non-self-similar black hole in a 
flat Friedmann background. In particular, the black hole 
singularity is connected to the big bang singularity and 
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z s 



z s 



z=z s 




z=z 2 



z=oo 

FIG. 5: The form of S{z) and V{z) for the asymptotically 
Friedmann solutions which contain a black hole with an event 
horizon, together with the conformal diagram. There is a fan 
a.t z — and a splash a,t z = z^h, and the causal structure 
is sSfTpSs. 



necessarily naked for a while. For E < E^, the mini- 
mum of V is above 1, as illustrated in Fig. ^, and the 
entire Universe is inside the black hole. (However, there 
is always an apparent horizon since CC show that the 
minimum of M is necessarily below 1/2.) In this case, 
the conformal diagram has two spacelike singularities and 
there are no null infinities. This is because the integral 
over (||), for both signs, converges as z ^ oo and z — > ^s, 
so any null geodesic covers only a finite range of In r be- 
tween the big bang singularity and the black hole/big 
crunch singularity. This contrasts with the conformal di- 
agram for a recoUapsing Friedmann universe, which has 
two separated spacelike singularities. 

Note that when the minimum of V exactly equals 1, the 
conformal diagram in Fig. |5| loses the part in which the 
homothetic lines are timelike and so there is no longer 
a naked singularity. In terms of causal structure, the 
"fTp" part of Fig. || is replaced by a "semi- fan" and a 
"semi-splash" . However, the conformal diagram is still 
very different from the one shown in Fig. ^ so there is a 
qualitative change as the minimum of V passes through 
1. In the cases discussed below, similar considerations 
apply whenever V has an extremum at ±1, but we will 
not comment on this explicitly. 



B. Asymptotically quasi-static solutions 

As discussed by CC, there is exactly one self-similar static 
solution for each value of a and a 1-parameter family of 
solutions that are asymptotically static (in the sense that 
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the radial 3- velocity Vr tends to zero as \z\ oo). There 
is also a 2-parameter family of solutions that are asymp- 
totically "quasi-static" (in the sense that Vr is finite but 
not necessarily zero). However, it should be emphasized 
that the solutions in this class may only be close to the 
static solution for part of their evolution. One of the two 
parameters can be taken to be the asymptotic energy E, 
while the other (denoted by D) gives the value of z at 
the big bang or big crunch singularity. Note that D, un- 
like E and Vr, cannot be expressed explicitly in terms 
of the asymptotic parameters, so it is sometimes more 
convenient to use Vr as the second parameter. 

The key feature of these asymptotically quasi-static solu- 
tions is that they span both negative and positive values 
of z and necessarily pass from z — ^oo to z = +00, 
whereas the asymptotically Friedmann solutions are con- 
fined to 2; > or z < and are symmetric. This is 
because the big bang occurs at z — —1/D (correspond- 
ing to a negative value of t) in these solutions, so the limit 
|z| — > 00 has no particular physical significance. These 
solutions can be interpreted as inhomogeneous cosmolog- 
ical models with an advanced big bang. Equivalently, for 
the time-reversed solutions, there is a big crunch singular- 
ity at z = +1/D. There are two types of asymptotically 
quasi-static solutions. Those with E exceeding some neg- 
ative critical value Ecr\t{D) expand or collapse forever. 
Those with E < Ecrit{D) expand and then recoUapse. 

1) Ever-collapsing solutions. These start out from an 
infinitely dispersed state and describe the collapse of an 
inhomogeneous gas cloud to a singularity at z = +1/D 
(i.e., after t = 0). (The ever-expanding solutions are 
just the time reverse of these and will not be considered 
exphcitly) They are of type TS or TSTS and illustrated 
in Fig. |7| and Fig. |. They start with y = at z = 
and then, as z decreases, reach a sonic point where V = 
^/a. In this context, it should be noted that the second 



FIG. 7: The form of S{z) and V{z) for the ever-collapsing 
asymptotically quasi-static solutions which contain a naked 
singularity, together with the conformal diagram. The causal 
structure is rTfSfTpSs. 



parameter D has relatively little effect on the form of the 
solutions in the subsonic regime. Indeed, all solutions 
apart from the exactly static one must be asymptotic to 
the flat Friedmann solution at small \z\. In particular, 
the models can collapse from infinity (i.e. S" — > 00 as 
z —f or t —00) only if E is positive or lies in discrete 
bands if negative. The subsonic solution is then attached 
to a supersonic asymptotically quasi-static solution at 
the sonic line. The supersonic solutions pass through a 
Cauchy horizon at zi (where V = 1) before tending to the 
quasi-static form at z = —00 and jumping to z = +00. 
As z further decreases, V first reaches a maximum and 
then diverges to minus infinity when it encounters the 
singularity ed z — 1/D. The maximum will be above -1 
if E is less than some negative value E+{D) and there 
will then be two points Z2 and Z3 where V — —1. This 
case is illustrated in Fig. |^ and the conformal diagram 
is a combination of the one shown in Fig. ^ and Fig. |^. 
This shows that one necessarily has a naked singularity, 
as pointed out by Ori and Piran. The maximum of V will 
be below -1 for E > E+{D) and this case is illustrated 



in Fig. |8 
in Fig. 



The conformal diagram now resembles the one 
except that the singularity has z = 1/D rather 



than an infinite value of z. 

2) Expanding and recollapsing solutions. Solutions with 
E less than the critical value £'crit(^) expand from an 
initial singularity in the z < region and then recoUapse 
to another singularity in the z > region. They are of 
type S, STS or STSTS. The value of z is at the 

initial singularity but depends on both E and D at the 
final one. As z decreases from —1/D, V rises from —00, 
reaches a maximum below — ^/a and then tends to the 
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FIG. 8: The form of S{z) and V{z) for the ever-collapsing 
asymptotically quasi-static solutions which do not contain a 
naked singularity. The conformal diagram is the same as in 
Fig. y, with structure sSfTr, although upside down. 
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FIG. 9: The form oi S[z) and V[z) for the recoUapsing asymp- 
totically quasi-static solutions which contain two singularities, 
one of which is naked. This shows the case Knax > — 1 and 
Vmin > 1. Both this case and the case with Knax < — 1 and 
Knin < 1 have the conformal diagram of Fig. ^ except that 
the big bang singularity is at jz = —1/D instead of 2: = oo. 



quasi-static form as z ^ —00. The solution then jumps 
to z = -l-oo and enters the z > regime. As z continues 
to decrease, V decreases to a minimum above y/a and 
then tends to -l-oo at the value of z corresponding to the 
recoUapse singularity. We now have four possible situa- 
tions, depending on whether the minimum of V is above 
or below +1 and whether the maximum of V is above or 
below —1. If Knax > — 1 and Knin > 1, one necessarily 
has a black-hole event horizon at zi and a cosmological 
particle horizon at Z2 in the z < regime. This occurs 
if E exceeds some negative value E^,{D) and the mini- 
mum of V will reach ^/a when E reaches i?crit(-D), cor- 
responding to the last recoUapsing solution. In this case, 
the solution is of the STS type and illustrated in Fig. ||. 
The conformal diagram just resembles that in Fig. |^, ex- 
cept that the initial singularity is at z = 1/_D rather 
than z = 00, and this shows that the black hole singu- 
larity is necessarily naked. The case with Kmax < — 1 
and Knin < 1 is equivalent to this and not shown ex- 
plicitly. The case with Knax > — 1 and Knin < 1 is of 
the STSTS type and illustrated in Fig. |l^. There are 
now four points (zi, Z2, Z3, Z4) at which K = 1 and so the 
conformal diagram is more complicated, though the sin- 
gularity structure is unaltered. The case with Knax < ^ 1 
and Knin > 1 is illustrated in Fig. |ll|. There are now no 
points where V — 1 and so the conformal diagram is 
as in Fig. ^. In principle, there might be solutions with 
E < E„it{D) in which V has neither a maximum nor a 
minimum. However, such solutions would have two sonic 
points and parameter-counting suggests that it is unlikely 
that any solution having crossed the sonic surface once 
in a saddle could cross it again. 

The exact static solution is illustrated in Fig. |l^. The 
conformal diagram has the Minkowski form at large dis- 
tance but there is a timelike naked singularity at the 
origin and points zi where K = — 1 and Z2 where K = 1. 
Note that there are no solutions which are asymptotic to 
the static solution at z = apart from the static solution 
itself, so there are no perturbations of the exact static 
solution. The asymptotically quasi-static solutions with 
a regular centre are asymptotically Friedmann at z = 0. 




z=z 4 



z=-l/D 



FIG. 10: The form of S{z) and V{z) for the recoUapsing 
asymptotically quasi-static solutions which contain two naked 
singularities, together with the conformal diagram. The 
structure is sSfTpSfTpSs. 



The asymptotically Minkowski solutions with 



There are two families of asymptotically Minkowski so- 
lutions in the z > regime and both are physical only 
for a > ^. Equivalent (time reversed) solutions exist in 
the z < regime but we do not discuss these explicitly. 

A) The first family is described by two parameters and 
has K — > 1, — > cxD and /i — )■ at some finite value 
z = z.^,. In this limit, 



V IV ^ ~{ba-l)/{l~ a). 



(13) 



This family is illustrated in Fig. 13. Although the value of 
z is finite at z*, the Schwarzschild radial distance (i? = 
rS) is infinite, as is clearly seen in the S{z) diagram. 
These are the "explosive solutions" of Ori and Piran [|], 
in which the fluid velocity approaches the speed of light, 
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FIG. 11: The form of S{z) and V{z) for the recoUapsing 
asymptotically quasi-static solutions which contain no naked 
singularities. The conformal diagram is the same as in Fig. ^ 
with structure sSs. 



s 




FIG. 12: The form of S{z) and V{z) for the exact static 
solution, together with the conformal diagram. The structure 
is sTfSfTs. 



leaving null data for flat spacetime on an outgoing null 
cone. 

B) The second family is described by one parameter and 
has V Vt, >1,S'— >cxD and /i — > as z — > cx). The 
expression for is 



approaches the Milne model along certain time lines. 

As can be seen in Figs. ^ and both families con- 
tain solutions that can be connected either to the origin 
z = (via a sonic point) or to a singularity (for which 
5 — > at some finite value zg). The former are of type 
T or ST and correspond to singularity-free models which 
collapse and then bounce at some finite density into an- 
other expansion phase. The latter are of type S and 
collapse to a singularity without passing through a sonic 
point. We therefore have four situations. The bounc- 
ing 2-parameter case is illustrated in the left conformal 
diagram of Fig. The conformal structure is very sim- 
ilar to the Minkowski case and can indeed be extended 
to a solution whose conformal structure is exactly like 
Minkowski if one adds another Minkowski patch to the 
top-left corner (shown dotted). There is a point singu- 
larity at r = t = 0. Since /ir^ = at z = z*, there is 
a singularity only where r — 0. The conformal diagram 
for the collapsing 2-parameter case is shown on the right 
of Fig. |l3|. This could also be attached to a Minkowski 
region (shown dotted), in which case it is similar to the 
conformal diagram in Fig. ^ It might also in principle 
be patched on to the top corner of the left conformal di- 
agram. The 1-parameter case is illustrated in Fig. 
Although V = Vit > 1 is finite as z — > oo, the surface 
z = cxD is not actually spacelike but consists of two null 
branches. As shown in appendix ^ one of them is an 
outgoing null surface at finite distance - as in the type 
A solutions, it is possible to glue it onto Minkowski here 
- and the other branch is /+. 

Note that there is a 1-parameter family of singular solu- 
tions for each value of zg. Most of these will be asymp- 
totic to either a type-A asymptotically Minkowski solu- 
tion or a quasi-static solution. However, for sufficiently 
large values of zs, there will also be one type-B asymptot- 
ically Minkowski solution and one asymptotically Fried- 
mann solution. The limiting value of zs is the same in 
each case, reflecting the fact that the conditions at large 
values of z have little influence on what happens near the 
singularity. 



Y{a + 1) + ^ya{a^ - a"^ + 3a + 1) 
1 — a 



(14) 



This family is illustrated in Fig. |lj. Note that the T4 
decreases as a decreases and reaches 1 when a — 1/5. 

In both cases, the limit S oo can be regarded as cor- 
responding to an infinitely dispersed state, analogous to 
the late stage of an open Friedmann model. Also both 
families of solutions have m/R in the asymptotic 
limit. The mass m itself tends to a finite value in case 
A but to zero in the case B. It is important to realize 
that these solutions are not asymptotically flat in the 
usual sense. Rather, they are perfect-fluid spacetimes for 
which the Minkowski geometry is obtained asymptoti- 
cally along certain coordinate lines. This is reminiscent 
of the open Friedmann solution, which asymptotically 



D. Asymptotic Kantowski-Sachs solutions 



The final class of models is associated with the 
Kantowski-Sachs solution. For each a there is a unique 
self-similar Kantowski-Sachs solution and there also ex- 
ists a 1-parameter family of solutions asymptotic to this 
at both large and small values of |z|. Solutions with 
— 1/3 < a < 1 are probably unphysical because the mass 
is negative and they are also tachyonic for < a < 1. So- 
lutions with —1 < a < —1/3 avoid these unsatisfactory 
features. Although such equations of state violate the 
strong energy condition, they could could well arise in 
the early Universe due to inflation or particle production 
effects. 
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FIG. 13: The form of S{z) and V{z) for the 2-parameter 
bouncing and collapsing asymptotically Minkowski solutions, 
together with the conformal diagram. There are two solu- 
tions which are geodesically incomplete. Each of them can 
separately be continued as Minkowski, and this is shown by 
the thin dashed line. One fluid world line is shown as a thicker 
dcished line. The two solutions can also be glued together at 
z = Zs- Note that in this case 2 = z* is a fan on one side and 
a splash on the other, because V{z) has a minimum eit V = 1 
there, rather than crossing V = 1. The causal structure of 
the incomplete spacetime on the left is rTfSfTf. The causal 
structure of the incomplete spacetime on the right is pSs. 




FIG. 14: The form of S{z) and V{z) for the 1-parameter 
bouncing and collapsing asymptotically Minkowski solutions. 
The conformal diagram and the causal structure are the same 
as in Fig. O, with z, replaced by 2 = 00. 



V. CRITICAL SOLUTIONS AND SHOCKS 

Critical phenomena in gravitational collapse arise during 
the evolution of asymptotically flat regular initial data 
that are close to the threshold between black hole forma- 
tion and dispersion. Such near-critical data approach a 
universal intermediate attractor before they either form 
a black hole or disperse. In dynamical systems terms, 
this "critical solution" is an attractor of codimension one, 
whose stable manifold separates the phase space into so- 
lutions that form a black hole and those that disperse. 
In terms of perturbation theory, it must have exactly one 
unstable linear perturbation mode, such that adding any 
amount of the unstable mode with one sign eventually 
drives the solution to collapse, and adding any amount 
with the opposite sign drives it to dispersion. 

The property of being at the boundary of collapse can be 
realized in two ways: either the solution is static, corre- 
sponding to an unstable star, or it is self-similar, corre- 
sponding to scale-invariant collapse to a point. In "type 
11" critical phenomena the critical solution is self-similar 
(either CSS or DSS), and the black hole mass mbh as a 
function of the initial data exhibits a power-law scaling 
of the form 



In the asymptotically KS solutions, S (x and V is also 
a decreasing power of z for —1 < a < —1/3. Thus both 
S and V decrease to as z ^ 00, as indicated in Fig. |l^. 
This means that there is a timelike singularity at t = 
and there is no spacelike infinity since R depends only 
on t. The conformal diagram is therefore as indicated in 
Fig. |l5|. The KS solution can be obtained from the static 
one by changing the value of a and interchanging r and 
t] this explains why the conformal diagram of the KS 
solution resembles one half of that of the static solution 



rrihh ^ [p-p*y ■ (15) 

Here p is any one parameter of the initial data that is be- 
ing varied while all other parameters remain fixed, and 
is its critical value, such that data with p > p^, form 
a black hole, while data with p < p^ disperse. The crit- 
ical exponent 7 is independent of the initial data, and 
depends only on the symmetry of the system under con- 
sideration and type of matter. It can be related to the 
growth rate of the unstable mode of the critical solution 
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Type II critical phenomena were first described by Chop- 
tuik [|ll| for the sphericaUy symmetric massless scalar 
field. Abrahams and Evans Q described type II criti- 
cal phenomena for axisymmetric vacuum gravity. In both 
cases the critical solution is DSS. The critical collapse of a 
spherically symmetric perfect fluid, the system which in- 
terests us here, was first investigated by Evans and Cole- 
man 1 13 and is CSS. Many other systems were investi- 
gated subsequently and it seems that critical phenomena 
are generic. For a recent review, see Q . 

Critical phenomena arise from any smooth initial data, 
in particular from analytic initial data. The time evolu- 
tion of these data is analytic until the singularity occurs. 
(In fluids, it is possible that a shock could form first, but 
empirically this is not the case.) The critical solution, 
if it is to approximate such evolutions, must share this 
analyticity property. A spherically symmetric CSS fluid 
critical solution in particular must be analytic at the cen- 
ter to the past of the singularity and at the sonic point to 
the past of the singularity. These requirements deflne a 
boundary value problem between the regular center and 
the first sonic point. This problem admits a countable 
family of solutions which are parameterized by the num- 
ber of zeros of the radial velocity Vr. The critical solution 
(with one unstable mode) turns out to be the one with a 
single zero ||l^ , with the fluid ingoing in a central region 
and outgoing in the exterior. These solutions were first 
obtained numerically in jl^ for a = 1/3, in ||l^, |l5j for 
< a < 0.89, and in [[6| for < a < 1. 

Type II critical solutions have renewed interest in the 
study of self-similar solutions in general relativity. We 
have seen that every self-similar spacetime has a strong 
curvature singularity s = (unless it is flat). In critical 
solutions this singularity is naked, in the sense that light 
rays can escape from regions of arbitrarily high curvature 
with a flnite redshift. In the collapse of near-critical ini- 
tial data, the solution is well-approximated by the critical 
solution down to a lower cut-off scale of s ~ |p— |. Con- 
sequently the maximal value of the Kretschmann scalar 
that can be seen is \p — p*\^'^- In the limit in which the 
parameter p is tuned to its critical value p*, the naked 
singularity of the critical solution is seen. Therefore the 
naked singularity of the critical solution is produced from 
a set of codimension one in the space of smooth and 
asymptotically flat data. 

Because the critical solution is determined by analyticity 
at the center in the past of the singularity and the past 
sound cone of the singularity, its continuation from the 
past sound cone is unique up to the future light cone of 
the singularity, which acts as a Cauchy horizon. We shall 
briefly discuss possible continuations beyond the Cauchy 
horizon, but this is likely to be irrelevant to gravitational 
collapse because quantum gravity effects are expected to 
become important at very high curvature. 

The nature of the critical solution in spherical perfect 
fluid collapse with the equation of state p = aji has been 



discussed by CCGNU @: For 0.28 < a < 1 the criti- 
cal solution is of the asymptotically Minkowski A type, 
and for the special value a ~ 0.28 it is asymptotically 
Minkowski type B. In both cases the spacetime structure 
is given by the left conformal diagram in Fig. 

The simplest continuation beyond the Cauchy horizon at 
z = z^, would be a wedge of flat empty spacetime, as 
indicated by the thin dashed line in the left diagram in 
Fig. ^ In this case, the CSS singularity is a single point. 
The Ricci scalar at the center diverges as t^^ for t < 0, 
and vanishes for < > 0. 

Another possible continuation would be the right con- 
formal diagram in Fig. |l^. From the general arguments 
in Section III we see that such continuations form a 2- 



parameter family. In this case the singularity has a null 
branch which is completely covered by a spacelike branch. 
Note that this continuation is filled by fiuid particles that 
emerge from the null branch of the singularity and end at 
the spacelike branch. It could therefore be described as 
a black hole expanding at the speed of light or as a baby 
universe. A distant observer would see the t^^ divergence 
of the curvature at the center, but would be swallowed 
by the expanding black hole at i = 0. 

For < a < 0.28, the critical solution is asymptoti- 
cally quasi-static, but its causal structure differs from the 
asymptotically quasi-static solutions that we have dis- 
cussed already. The reason is that it crosses the sonic 
surface twice. The structure of the critical solutions with 
< a < 0.28 up to the second sonic point is TST and 
shown in Fig. We call the second sonic point the sonic 
horizon, in analogy to the Cauchy horizon of the singu- 
larity. CCGNU find that at the sonic horizon the critical 
solution reaches the sonic surface away from the sonic 
line. As discussed in Appendix this means that it 
cannot be continued as a C*^ solution, so the density gra- 
dient and velocity gradient become infinite at the sonic 
horizon. 

This suggests that before the solution curve reaches the 
sonic surface a shock occurs. It must be to the future of 
the Cauchy horizon and to the past of the z-line where 
the CSS ansatz breaks down. Therefore its location curve 
r — fsh(i) is bounded by Zch < rsh{t)/t < Zgp. The most 
natural assumption compatible with this is that the shock 
location is also self-similar, or rsh_{t)/t = Zsh — const. 
Such solutions have been investigated by Cahill and Taub 
|]l7| , and we need only summarize their results here with 
our notation. 

Given a CSS solution with a shock at zq, the solution on 
one side uniquely determines the solution on the other 
side. The Israel matching conditions reduce in spherical 
symmetry to continuity of the Hawking mass m and the 
area radius i?, so it follows that M — m/ R is continuous. 
We can also use the scale invariance of CSS solutions to 
impose the gauge condition that the coordinate r be con- 
tinuous. Continuity of R is then equivalent to continuity 
of S. The Rankine-Hugoniot conditions for the fluid give 
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FIG. 16: The form of S{z) and V{z) and the conformal dia- 
gram for the critical solutions with < a < 0.28. The space- 
time cannot be continued through the second sonic point as a 
CSS and continuous solution, but it can be continued through 
a shock. The two dashed lines are the sonic points, marked 
with dots in the S and V plots. 



V and 



SttG 



(16) 



on one side in terms of their values on the other side (Eqs. 
(6.23-24) of [|7j). The complete matching conditions are 



S+ 
M+ 

w+ 



(y2 



(17) 
(18) 

(19) 
(20) 



Note that \V\ < ^Ja on one side and V > ^Ja on the 
other. One also requires |y| > a on both sides for W to 
be positive. Therefore V- is restricted by the condition 
a < \y~\ < \/a. We thus have a 1-parameter family of 
possible continuations parameterized by V- in this range, 
so the existence of the shock adds one free parameter, 
namely zq, just like a nodal sonic point. In the strong 



shock limit V_ 



we obtain W+ and Vr — > 1, 



but the energy density measured by a timelike observer, 
/i+/(l — y^), remains finite in this limit. In the limit 
vl ^ a the shock disappears, with W+ ~ W- and V+ = 
V-. 



VI. CONCLUSIONS 



In this paper we have established the conformal struc- 
ture of all possible spherically symmetric CSS solutions, 
independently of the Einstein equations. We have then 
listed all possible conformal structures if the matter is 
a perfect fluid with equation of state p — a/i. This in- 
cludes the solutions associated with critical phenomena 
in gravitational collapse. 

Our kinematical discussion is based on the ideas of Nolan 
Q and some of our diagrams were first drawn by Ori and 
Piran ||l|. However, Nolan's kinematical discussion is re- 
stricted to spacetimes compatible with ingoing null dust 
matter, and his use of the Einstein equations is limited to 
imposing a positive energy condition. Ori and Piran only 
consider spacetimes with a regular center in the past. 

The main limitation of our discussion is that we have ex- 
cluded solutions with shocks. Such solutions may have 
more free parameters than diffcrcntiable ones and they 
may even admit additional conformal diagrams. Also we 
have only considered CSS continuations of the critical so- 
lutions beyond their Cauchy horizon but this is probably 
not a serious restriction because the continuation must 
be asymptotically CSS at small scales. 

By drawing the conformal diagrams we have been able 
to elucidate the nature of some previously obscure CSS 
solutions. In particular, the conformal structure of the 
asymptotically Minkowski solutions of type A and type 
B are the same, and the CSS solutions that describe a 
black hole forming in an open universe have a null sin- 
gularity, which can either be completely covered by the 
black hole singularity (see Fig. ^on the right) or covered 
only partly. 

An interesting result is that, while the singularity may 
have past branches from which fluid worldlines emerge 
and future branches where they end, these branches must 
be connected by a null branch. All parts of the singularity 
are therefore topologically connected. The recoUapsing 
models, for example, do not have the conformal diagram 
of the closed Friedmann universe, with disjoint big bang 
and big crunch singularities, but one of the conformal 
diagrams indicated in Fig. ^, Fig. ^ and Fig. |l^. 

Reviewing earlier work ||^, we have placed the critical 
solutions which arise in gravitational collapse within our 
general classification. These solutions are analytic and 
unique up to the Cauchy horizon of the CSS singular- 
ity. Intuitively one might expect the continuation of the 
critical solution to be non-unique because it depends on 
information coming out of the naked singularity. As in- 
formation in spherically symmetric fluid solutions travels 
only at the speed of sound (because there are no grav- 
itational waves), one would expect the spacetime to be 
non-unique beyond the sonic horizon of the singularity. 
However, contrary to this intuition, we have shown that 
the critical solution is non-unique even before the sonic 
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horizon, at the Cauchy horizon. 
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^2, while the flat Friedmann solution has a sonic point 
above Z3. 

The two nodal segments of the sonic line are attractors, 
and are crossed by generic 2-parameter families of solu- 
tions. The saddle segment of the line is crossed only by a 
1-parameter family of solutions. Generic solutions there- 
fore miss the saddle segment of the sonic line and reach 
the sonic surface away from the sonic line. Expanding 
around a point in the sonic surface but off the sonic line, 
where — a but .g 7^ 0, we find to leading order that 



APPENDIX A: THE SONIC SURFACE 



Bickncll and Henriksen Q have expanded a generic so- 
lution around a point in the sonic line in terms of a pa- 
rameter u along the solution curve (M, S, z^rf). To linear 
order they find 



{6M, 5S, dz, S'i])(u) = aie 



a2e 



\iu-t 



(Al) 



where the eigenvalues Ai 2 and eigenvectors Vi 2 are 
known, and a 1,2 are free constants. By parameter count- 
ing, this linearized solution is generic. There are also 
two other eigenvectors with zero eigenvalues in the 4- 
dimensional {M,S,z,r]) space. One takes the solution 
out of the constraint surface c = 0, and the other takes 
it along the sonic line. Their amplitudes can consistently 
be set to zero. 

The nature of the solution depends crucially on Ai^2, 
which are the solutions of a real quadratic equation. If 
they form a complex conjugate pair the point on the 
sonic line is called a focus. The formal solution (Al) 
is then complex and spirals around the sonic point in 
phase space. No solution crosses the sonic surface at this 
point. 

If A1.2 are real with the same sign, the point on the sonic 
line is called a node. We may assume without loss of 
generality that Ai < A2 < 0. Then the sonic line is 
crossed as u ^ 00. Furthermore, as u — > 00, e^^" can 
be neglected compared to e'^^", and the implicit functions 
M{z), S(z), r]{z) are C^. Therefore a 1-parameter family 
of solutions crosses each point of the sonic line. All 
these solutions have the same tangent at the sonic point. 
The solution with ai = from this family is analytic. 
There is also an isolated analytic solution with 0:2 = 0. 

If Ai^2 have opposite signs, the point on the sonic line is 
called a saddle. Without loss of generality, we assume 
Al < < A2. Then there are only two isolated solutions 
that cross the sonic surface: the solution with ai = as 
u —00 and the solution with a2 = as u — * 00. Both 
solutions are analytic. 

The sonic line can be parameterized by z. It then turns 
out that the points < z < zi are saddles, zi < z < Z2 
are nodes, Z2 < z < z^ are focal points, and Z3 < z < 00 
are nodes again, where zi_2.3 are known functions of a 
PI . The static solution has a sonic point between zi and 
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±J2g{^] 6z, 



drj 



(A2) 
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(A4) 
(A5) 



This point on the sonic surface is reached by two solution 
curves from one side, but there is no solution curve that 
reaches it from the other side. Note also that, because 
R — rS and SttG^ — r^^r], the comoving physical density 
gradient d^/dR and velocity gradient dV/dR diverge. 



APPENDIX B: ASYMPTOTICALLY 
MINKOWSKI TYPE B 



From Eq. (3.43) of we see that as z ^ 00, the space- 
time metric is approximately of the form 



ds^ 
where 



-A^z^^dt^ + B^z^^'dr^ 



syz^- 



dn^ 



(Bl) 



B/A = V,, a- 1 = &= 1/(^2-1), c=l/(K-l). 

^ . . ^^^^ 

where 14 is given by Eq. (14). This (approximate) metric 

can be written as 



ds^ 



-C du dv 



where C > is a known constant, and 



u = -rz-^ni+v.) 



(B3) 



(B4) 



are negative and positive, respectively. Contrary to the 
claim in this metric is not Minkowski. From the form 
of the metric, it is clear that u and v are affine parameters 
of null geodesies, z = 00 has two branches: u = 0, 
which is therefore at finite affine parameter distance, and 
V = 00, which is at infinite affine parameter distance. 
The former is a regular surface and the spacetime can 
be continued there. The latter is /"'". Fluid worldlines 
r = const end up at (m = 0, u = 00). Note that this point 
is not ig: observers that resist the outward fluid pressure 



14 



can cross the surface u = 0. It should be stressed that R = Sqv. However, non-leading terms actually show that 
metrics ( |Bl[ ) and ( |B3| ) are only approximate. Were they m/R — > in the limit z oo. 
exact, the mass function would be m/R = 1/2 because 
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